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We establish a new covariant phenomenological model, which describes an inhuence of pseu¬ 
doscalar (axion) held on spins of test massive particles. The model includes general relativistic 
equations of particle motion and spin evolution in background pseudoscalar (axion), electromag¬ 
netic and gravitational helds. It describes both the direct spin-axion coupling of the gradient type 
and indirect spin-axion interaction mediated by electromagnetic helds. Special attention is paid to 
the direct spin-axion coupling caused by the gradient of the pseudoscalar (axion) held. We show 
that it describes a spin precession, when the pseudoscalar (axion) held is inhomogeneous and/or 
non-stationary. Applications of the model, which correspond to the three types of four-vectors 
attributed to the gradient of the pseudoscalar (axion) held (time-like, space-like, and null), are 
considered in detail. These are the spin precessions induced by relic cosmological axions, axions 
distributed around spherically symmetric static objects, and axions in a gravitational wave held, 
respectively. We discuss features of the obtained exact solutions and some general properties of the 
axionically induced spin rotation. 

PACS numbers: 13.88.-|-e, 14.80.Va, 95.30.Sf 


I. INTRODUCTION 

The terms angular moment and spin are consid¬ 
ered in theories of gravity in two main aspects. First, we 
face with them, when we study the gravitational field 
formed by a single rotatin g e xtended body (e.g., black 
holes, neutron stars, etc. li), or a system of bod¬ 
ies rotating around an attracting center (e.g., stars in 
spiral galaxies Si). Second, we deal with the polar¬ 
ization or spin, when we investigate the gravitationally 
induced dynamics of test particles, which possess vec¬ 
tor (or pseudo-vector, respectively) degrees of freedom. 
This sector of investigations is usually connected with 
experiments in gravitational physics (see, e.g., Si and 
references therein). 

According to the standard terminology the ro¬ 
tating extended body or point-like particle can 
be described by an anti-symmetric tensor of to¬ 
tal moment Sik, which can be decomposed as 
Sik=5^'^LmUn-(tikmnS^U'^, using the unit velocity 
four-vector [/" of the body. The quantity L™ can be in¬ 
dicated as a four-vector of the orbital moment, and the 
(pseudo)vector S'™ plays the role of a classical spin four- 
vector. The term is the four-indices Kronecker ten¬ 
sor and eikmn is the Levi-Civita (pseudo)tensor. Evolu¬ 
tionary equations for the total moment Sik and/or for 
its constituents L™ and S™ have been investigated in 
various contexts by many authors. The list of obtained 
results is very long, and we would like to attract the 
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attention to the following two ones only. 

In 1959 Bargmann, Michel and Telegdi obtained 
the covariant equations for classical spin particles with 
anomalous magnetic moment Q- This model (BMT- 
model, for short) describes a spin evolution in the 
framework of Special Relativity as a generalization of 
the non-relativistic theory elaborated by Thomas 


; orv e 

Frenkel [9| and Bloch [T^. Being completed by the 
spin-curvature coupling terms, introduced earlier by 
Mathisson m and Papapetrou [l^ , this model became 
a starting point for numerous investigations of s pin p ar- 
ticle dynamics in General Relativity (see, e.g., [l3l - l^ 
and references therein). When it was necessary to in¬ 
clude electromagnetic interactions, these investigations 
used the Faraday-Maxwell electrodynamics. Scientific 
events of two last decades attracted the attention to 
its generalization called Axion Electrodynamics. In the 
light of the hypothesis about axionic nature of the dark 
matter this theory seems to be more appropriate for de¬ 
scribing the cosmic electrodynamics than the standard 
Faraday-Maxwell theory. 

The story, how axions were associated with dark 
matter particles, is well-known. The dark matter, a cos¬ 
mic substance, which neither emits nor scatters the elec¬ 
tromagnetic radiation, is assumed to accumulate about 
23% of the Universe energy. The mass density distri¬ 
bution of the dark matter is presented in the astro- 
physical catalogues as a result of observations and the¬ 
oretical simulations (see, e.g., [IMl for details, review 
and references). The origin of the dark matter is not 
yet established. One of the most attractive hypothesis 
links the dark matter with massive pseudo-Goldstone 
bosons. The particles of this type were postulated in 
1977 by Peccei and Quinn in order to solve the 
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problem of strong CP-invariance, and were introduced 
into the high-energy physics as new light bosons by 
Weinberg 27| and Wilczek in 1978. Later these 
pseudo-bosons were indicated as axions and now they 
are considered as the most appropriate candidate for 
dark matter particles (29l - [^ . 

The description of these axions in terms of the field 
theory is based on the introduction of a pseudoscalar 
field (p, which is assumed to interact with electromag¬ 
netic and SU (2) - symmetric gauge fields. The theory of 
interaction between electromagnetic and pseudoscalar 
fields was elaborated by Ni in 1977 [s^. Due to the 
works of Sikivie (see, e.g., [s^) we recognize now this 
sector of science as axion electrodynamics] a number 
of effects such as the axionically induced optical activ¬ 
ity [ 13 , Ell magneto-electric cross-effects [H, etc., 
were predicted in the framework of this theory. 

Thus, on the one hand, the axions are associated in 
our mind with the cosmic dark matter, which regulates 
key gravitational processes in our Universe, for instance, 
provides specific (flat) distribution of star velocities in 
the spiral galaxies (see, e.g., H). On the other hand, 
the axions produce the photon polarization rotation, in 
other words, the axionic dark matter forms an electro¬ 
dynamically active chiral cosmic medium. Clearly, we 
could try to find fingerprints of the axionic dark mat¬ 
ter in the electromagnetic signals associated with astro- 
physical and geophysical phenomena [45l - l43| . 

Returning to the problem of influence of the (ax¬ 
ionic) dark matter on the orbital moment, polarization 
and spin, we can say the following. First, definitely, the 
dark matter influences the orbital moment of star gas 
in spiral galaxies, providing the observed flat profiles 
of the velocity distribution Q. Second, the axionic 
dark matter is predicted to rotate the photon polariza¬ 
tion [iO]. Third, due to the first and second arguments, 
one can expect that the axionic dark matter should also 
influence the particle spin. 

In fact, the theoretical discussions about a new 
branch in the axion physics, i.e., the theory of spin- 
axion coupling, were opened in 80s of the last century; 
it was a discussion concerning the axionically induced 
spin dependent long range forces [HO]. The long range 
spin dependent forces were in a focus of experimental 
investigations (see, e.g., [5l|), and this circumstance 
gave the impetus for development of new experimen¬ 
tal programs devoted to detection of axions, based on 
their interactions with nuclear spinSjClectric and mag¬ 
netic moments, etc. (see, e.g., (53 - 1601 1 and references 
therein). This branch of axion physics was indicated as 
spin-axion coupling by analogy with the axion-photon 
and the axion-gluon couplings. Nevertheless the term 
spin-axion coupling is rather wide and two theoreti¬ 
cal approaches have to be distinguished. Experimental 
works in the terrestrial laboratories deal with bonded 
spin particles in the specific material media, and these 
bonded particles are, in average, at rest with respect 
to these media. For these purposes the non-relativistic 
formalism of condensed matter physics is adequate and 


describes correctly the models of spin dynamics. When 
one deals with free spin particles, which move with high 
velocity (e.g., in storage rings) or in a vicinity of sources 
of strong gravitational and axion fields (e.g., near axion 
stars), the covariant formalism of high energy physics 
is necessary for description of the model. 


From the mathematical point of view, our goal 
is to develop the covariant version of the theory of 
spin-axion coupling, and we consider extended dy¬ 
namic equations for relativistic spin-particle. The 
corresponding extensions include axionic modifications 
of the Bagrmann-Michel-Telegdy force-like term, and 
gradient-type modifications of a rotatory term. From 
the physical point of view, adding the third entity, the 
axion, to the interacting pair of photon and charged 
spin particle, one can expect that three channels of in¬ 
teraction will be activated: first, the already known 
direct axion-photon coupling; second, the direct spin- 
axion coupling; third, the indirect spin-axion coupling 
mediated by the axion-photon interaction. In this pa¬ 
per we intend to discuss in detail two last channels of 
the photon-axion-spin interactions. 


The paper is organized as follows. In Section HIl we 
analyze some details of the well-known model of axion- 
photon coupling to motivate our new model of spin- 
axion coupling. In particular, we derive an equation of 
the photon polarization precession induced by a gradi¬ 
ent of the pseudoscalar (axion) field. In Section imi we 
generalize the results of Section [IT] for massive particles 
and reconstruct phenomenologically equations of spin- 
particle dynamics in the axion environment. We divide 
the obtained terms into three types: first, the terms 
indicated as axionic generalizations of the Bargmann- 
Michel-Telegdy terms; second, the terms describing di¬ 
rect spin-axion coupling; third, the nonminimal terms. 
In Section|lV]we consider four applications of the model 
with the direct spin-axion coupling. In Subsection IIV Al 
we obtain an exact solution to the dynamic equations 
for the massive spin particle, which moves straightfor¬ 
wardly and interacts with the relic cosmic axions. In 
Subsection IIV Bl we apply the formalism to the model 
of relativistic charged spin particle motion in a storage 
ring with magnetic field. In Subsection IIV Cl we discuss 
the spin-axion coupling in the static spherically sym¬ 
metric gravitational background. We studied in detail 
two cases, describing radial and circular particle mo¬ 
tion, respectively. The exact solutions for the spin pre¬ 
cession in the field of pp-wave gravitational field pro¬ 
vided by axions are presented in Subsection IIV PI In 
Section El we discuss the common details of the axion¬ 
ically induced spin precession, which were revealed in 
all three submodels. 
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II. CLASSICAL ANALOGY: PHOTON 
POLARIZATION ROTATION INDUCED BY 
THE AXION FIELD 

In 1977 Ni [s^ obtained equations for the electro¬ 
magnetic field coupled to a pseudoscalar field (j). The 
initial form of the currentless electrodynamic equations 
was the following: 


as fc™ = h , where r is an affine parameter along 
the line, and h is the Planck constant. 

The differential consequence of the eikonal equa¬ 
tion kmk"^ = 0 is k^Vnkm = 0. Since the quantity fc"* 
is the gradient four-vector, the relation Vrnkn = '^nkm 
is valid, so we obtain additionally that 

= ( 8 ) 


= 0 . ( 1 ) 

This equation is derived by variation with respect to 
the potential four-vector Ai of the electromagnetic part 
5 '(em) of the total action functional, where 



prnn , pr. 


■ ( 2 ) 


Here F^n is the Maxwell tensor, and F*’^^ = 
is its dual tensor. As usual, Levi- 

Civita tensor, based on the completely anti-symmetric 
Levi-Civita symbol with The Maxwell 

tensor Fmn is connected with the electromagnetic po¬ 
tential Ak in the regular way: 


Tmn — ^mAji ^jiAjji , (3) 

SO, that the well-known relationship 

= 0 (4) 

converts into the identity. Due to (H]) the basic equation 
o simplifies as follows: 

ViF^^ =-F*^^Vi4). (5) 

In order to transform Eq. ([5]) to an one-photon form 
we use the standard procedure. First, we follow the 
approximation of geometrical optics, and represent the 
potential A^ as 


Am = amC*'*'. (6) 

Here is a rapidly varying phase, the so-called eikonal. 
Its gradient four-vector kj = V j 4' has a sense of a wave¬ 
length four-vector of the photon. The amplitude Um and 
its derivatives are assumed to be slowly varying func¬ 
tions, i.e., max|Vjam| min|fcjam|. In these terms 
the Lorentz condition 


For an illustration of main idea we assume that the 
wavelength vector A:™ is divergence-free, i.e., = 0 

(or equivalently, V™Vm'I' = 0). This can be realized 
in many models, e.g, in the Minkowski space-time, in 
the space-time with plane-wave symmetry. Then the 
hrst order approximation in Eq. ([5]) gives the following 
equation for the evolution of the amplitude : 

k’^Vma^ = Upkg . (9) 

The convolution of this equation with aj yields 

= 0 , ( 10 ) 

i.e., a^Qj = —a^ is constant along the photon world¬ 
line with the tangent vector kK This fact allows us to 
introduce the unit space-like polarization four-vector 
as follows: 


and to exclude the amplitude a from Eq. i 

Using the covariant differential operator 
D dx^„ 


Dr 


dr 


( 11 ) 


( 12 ) 


(13) 


one can rewrite Eqs. (|S]) and (nm as a pair of basic equa¬ 
tions of photon evolution: first, the dynamic equation 


Dt 


= 0 , 


and second, the equation of polarization rotation 

Dt 2 dr 


(14) 


(15) 


'^ruA^ = 0 (7) 

gives two relations. The first one, Umk"^ = 0, appears 
in the first order approximation and means that the po¬ 
larization is orthogonal to the wavelength four-vector. 
In the next order we obtain WmoF = 0. The leading 
order approximation in Eq. in the context of ® 
yields kmk^ = 0, the ordinary eikonal equation, imply¬ 
ing that the photons travel along null geodesics. The 
tangent vector to the photon world-line can be defined 


It is clear, that in terms of one-photon description three 
scalar quantities, k'^km, km^'^, and are the in¬ 

tegrals of motion for the set of equations (HI, (HU), 
i.e., they remain constant along the photon world-line, 
k'^km =const= 0, kmC^ =const= 0, = —1- 

There is an obvious analogy between equations of 
motion for the massless photons and massive fermions. 
Based on this analogy, we can replace the polarization 
four-vector y by the spin four-(pseudo)vector , in¬ 
troduce the time-like particle momentum four-vector 
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= mU^ = instead of the null wave four-vector 

(here and below we consider the system of units, in 
which c = 1). As a natural extension of equations dT41) 
and (|15l) we obtain the following evolutionary equations 
for the fermion particle: 


Dpi 

Dt 


= 0 , 


DSi 

Dt 


a 

2m 


J’- 


SP p^ 


(16) 


where a is some parameter introduced phenomenologi¬ 
cally. The presented analogy could explain the appear¬ 
ance of the rotational term in the right-hand side of the 
equation (HU). However, some fermions possess electric 
charges and thus interact with electromagnetic fields, 
and a phenomenological generalization of equations (jl4l) 
and m requires more sophisticated efforts. We will 
take into account these interactions in the next section. 


III. SPIN-PARTICLE DYNAMICS IN AN 
AXION ENVIRONMENT 


In the minimal theory the quantities T'' and Q'’ depend 
on the particle momentum spin S'*, Maxwell tensor 
Fmn and its dual as well as, on the pseudoscalar 

(axion) field (j) and its gradient four-vector In the 

nonminimally extended theory the quantities and 
t/* can include the Riemann tensor Ricci tensor 

Rik and Ricci scalar i?; also the covariant derivatives of 
the Maxwell tensor, S/sFmm and of the Riemann tensor 
^sRikran Can appear in the decomposition of F'" and 
C/®. Three equations (dll), (EOI) and m are satisfied 
identically, when 

= oj^'^Pk , = oj^'^Sk , ( 21 ) 

where w®*' is an arbitrary anti-symmetric tensor, = 
As a hrst step, we remind to the reader, what is 
this tensor for the most known example, the Bargmann- 
Michel-Telegdi model Q. 


A. Basic equations 


B. Bargmann-Michel-Telegdi model 


Let us consider the evolution of relativistic point 
particle with an electric charge and a spin four-vector. 
Dynamic equations for the particle momentum p® and 
for the spin four-vector S'* can be written as: 



PS* 

Dt 


(17) 


i.e., the rates of change of the given quantities are pre¬ 
determined by the corresponding force-like terms, F^ 
and t/*, respectively. Let us present, first, three general 
properties of the four-vectors J"® and t/*. 

(i). The mass of the particle, m, defined from the nor¬ 
malization law p*pi = w?(?, is assumed to be conserved 
quantity, providing the four-vector F'' to be orthogonal 
to the momentum: 


Dv^ 1 D 

= (18) 

(ii). Similarly, we assume that the scalar square 
of the space-like spin four-vector is constant, i.e., 
S^Si =const= —5^. Then using the second equation 
from dm we obtain that 


n c* ID 


1 D 


5^ 


0, (19) 


or in other words, that the force-like term C/* is orthog¬ 
onal to the spin four-vector. 

(in). Finally, we assume that the spin four-vector is 
orthogonal to the momentum four-vector, S^pi =const= 
0; then one obtains that 

^(p*.s,) = 0 ^ rs, + G^f = 0. (20) 


Let us consider the Bargmann-Michel-Telegdi 
model, for which the relevant force-like term F^ is the 
Lorentz force 



( 22 ) 


and G'" is of the form 


G^ 


e 

2m 


9F\S 


k 


(g-2) 


P^FkiS’^p^ 


(23) 


Here g denotes the so-called p-factor; the case g ^ 2 
means that the particle possesses an anomalous mag¬ 
netic moment. Clearly, these terms can be written in 
the form (I^Tll . when 


(jj 


ik 


= UJ 


ik _ ^ 

2m 



(g-2) 

2 




(24) 


where <5^^^ = 5(^5^ — is the four-indices Kronecker 
tensor. The important detail of this model is that the 
force F^ does not contain the spin four-vector, and the 
term C/® is linear in S^. There is a simple motivation 
of such model construction. In the quasi-classical ap¬ 
proach one uses the decomposition of microscopic equa¬ 
tions with respect to small quantity h, the Planck con¬ 
stant. Although the spin of particle enters the micro¬ 
scopic equations in the product h ■ S, the Planck con¬ 
stant is not involved into the classical dynamic equa¬ 
tions. Therefore one has to exclude the quantity h ■ S 
from (l24l) to provide that both left-hand and right-hand 
sides of Eqs. (II3 are of the same order in H, and the 
multiplier h can be eliminated. 
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C. Electrically charged spin particle in 
electromagnetic and axion fields 

1. Reconstruction of the force-type sources 

Keeping in mind the general relationships (I21|) we 
reconstruct the tensor using the following ansatz: 

A) The tensor does not contain the four-vector 5*. 

B) The tensor w** is up to the first order in the Maxwell 
tensor 

C) The tensor is linear in the pseudoscalar (axion) 
field 4>, or in its gradient four-vector 

D) The tensor is linear in the Riemann tensor. 

The first two requirements are the same as in the 
Bargmann-Michel-Telegdi model. The third point is a 
new detail, which appears as a natural extension of this 
model for the case of particles interacting with the pseu¬ 
doscalar field. The last point implies that the model can 
be minimal, when does not enter the tensor w®^, 

and nonminimal, when there are terms containing the 
Riemann tensor and its linear convolutions. Below we 
consider all the appropriate constructions, which satisfy 
the requirements A)-D) and can be added to the tensor 

(IMl)- 

Minimal (curvature free) terms linear in the axion 
field (j) itself can be represented as follows: 


procedure. First, we replace the Maxwell tensor F®^ 
in (IMl), dMD, and ()26|) with the tensor of nonminimal 
polarization-magnetization M®^ = 7?.®^'®®"Fmn, where 
the so-called nonminimal susceptibility tensor is 

introduced according to the rule 

^ikran ^ - 5®®®^'="®) -p ggR®'®®®®'® 

_ RiUgkm J^kng^m _ J^km ^ 3 ?) 

Here i?®®®®® is the Ricci tensor, R is the Ricci scalar, 
9i, q?., q.s are the nonminimal coupling constants (see, 
e.g., [Gll - I^ for details). Similarly, we replace the dual 
Maxwell tensor F*®* with the tensor A^®^: 

Al®'® = [ai *7^®'®®"’® -b aaF*®'®®®®’®] F„„ , (28) 

where *7^®*^®®®®® and 7^*®^®®®®® are the left-dual and right- 
dual tensors of nonminimal susceptibility, respectively. 
The corresponding coupling parameters uiab should be 
replaced with ojab- Concerning the last term in (1261) . 
which does not contain the Maxwell tensor, our strategy 
is to add the following nonminimal analogs: 

Knm) = {(h25 [di P, + 

+cc26p'e*'=^„7^®®®»4 . (29) 


ik 




• (25) 


In fact, can be obtained from with replace¬ 
ments Fmn and g ^ gA, where the coupling 

constant gA is an axionic analog of p-factor. Dimension¬ 
less parameter A is equal to one, if the Nature admits 
this coupling term, and A = 0, if it does not admit. 

Minimal terms linear in the gradient four-vector 
of the axion field can be written as follows: 


Here 72.’®®'^ = 72®®®®®'*'*(?„s. As a result of the described 
procedure, we deal with a large number of phenomeno¬ 
logically introduced coupling parameters. As a first 
step, below we restrict ourselves by the minimal model, 
i.e., by the case when ujab = 0; the nonminimal model 
contains a lot of specific details and will be discussed in 
a special paper. 


8. Dynamic equations 


ik 

W(2) 






goF 


*ik I tdG 2) 




(26) 




Here, in addition to the axionic analog of the p-factor 
gA, we introduced its gradient-type analog go- The 
constant /i in (IMl) plays the same role as the constant A 
in dMl). Other constants have no direct analogs, and we 
indicated them as a; 23 , W 24 , etc., where the first index 
2 is an indicator that the decomposition relates to the 
term W( 2 )- term in (IMl) is a unique element 

of the presented irreducible decomposition, which does 
not contain the Maxwell tensor F®®®®®. Note, that this 
term with a ;25 in front corresponds to the right-hand 
side of Eq. dHl) with the multiplier a. 

In order to classify nonminimal terms in the the¬ 
ory of spin-axion coupling we can use the following 


Following the representation of the total tensor 
cu®*^ = + w)*)) described above, we obtain 

the minimal dynamic equation in the form 

T)t}^ p 

^ - {F®'=pfc + F*®'=pfe [Ac/. + g (p'Vic/.)] } . (30) 

The first term in the right-hand side of this equation 
is the usual Lorentz force. The second term can be 
interpreted as an axionic analog of Lorentz force, in 
which the Maxwell tensor F®^ is replaced by its dual 
tensor F*®^, and the electric charge e is replaced by 
an effective pseudo-charge ecj). If A = 0 identically, 
the model is characterized by vanishing pseudo-charge. 
The coupling constant pi introduces a completely new 
term, which contains the gradient four-vector of an ef¬ 
fective pseudo-charge ecj). We faced with similar situ¬ 
ation in the 51/(2) and 5t/(3)-symmetric gauge mod¬ 
els, in which the isospin and color charges, respectively, 
were considered as functions, and their derivatives also 
entered the master equations [6^ . 
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The evolutionary equation for the spin with the 
redefinitions, mentioned above, takes the form 


DS^ 


Dt 


e 

2m 


gF^^Sk + 


+ 


eA 

2 m*' 


nm ^ 




gGF*^^Su + ,SV 




The first term in the right-hand side of this equation 
is the usual term attributed to the Bargmann-Michel- 
Telegdi model with anomalous magnetic moment; other 
terms can be indicated as its axionic analogs. Only one 
term with a ;25 in front does not contain the Maxwell 
tensor thus describing the direct spin-axion coupling. 


3. Our further strategy 


third example corresponds to the case = 0 , 

which can be realized in space-times with plane-wave 
symmetry (gravitational waves). In all three cases we 
consider the spin particle as a test one. It moves in a 
given space-time through the pseudoscalar field, which 
obeys the equation 

+ = (33) 

Here V((/)^) is the potential of the pseudoscalar (axion) 
field. This equation is derived from the axionic part of 
the total Lagrangian 

^(axion) = J ^ [V^cf^) - Vfc^VV] , (34) 

where the constant 4^0 is reciprocal to the coupling con¬ 
stant of the axion-photon interaction PA 77 , i-e., = 

PA 77 (see, e.g., 0). Since the direct spin-axion effect 
appears if and only if the gradient four-vector is 

non-vanishing, i.e., Vfc^ yf 0 , we assume that just the 
gravitational field produces the inhomogeneity or non- 
stationarity of the axionic field (j). 


We have established the phenomenological model, 
in which three channels of spin-axion coupling can be 
distinguished. The first channel is direct: it works even 
if there are no electromagnetic fields {Fmn=0), and the 
tidal (nonminimal) interactions are absent (72.*^"*" = 0 ). 
The second channel is indirect, the corresponding spin- 
axion coupling is mediated by the Maxwell field {F^n ^ 
0, 7?.®^™"=0). The third channel is also indirect, and it 
can be opened when the model is nonminimal ^ 

0). Below we put Fmn=^ in (15X1) and consider effects 
of a direct spin-axion interaction only. Next papers will 
be devoted to a systematic study of effects mediated 
by electromagnetic fields of various structures. In the 
future we also intend to add nonminimal couplings to 
the direct and indirect models of spin-axion coupling. 


IV. APPLICATIONS OF THE MODEL WITH 
DIRECT SPIN-AXION COUPLING 


We focus now on the direct spin-axion interactions, 
i.e., we assume that F^n = 0 and 7 ^®'="*" = 0 , and 
the gradient four-vector is non-vanishing due to 
the coupling of axions to the gravity field. In this case 
Eqs. o and (1311) obtain the simple form 



DS^ 

Dt 




(32) 


Below we study three examples. The first one relates to 
the cosmological context, and the gradient four-vector 
is time-like (i.e., > 0)- The second ex¬ 

ample relates to the case with the space-like gradient 
four-vector (i.e., < 0 ), which can be realized, 

e.g., in a spherically symmetric static space-time. The 


A. Spin coupling to relic dark matter axions 


Dark matter hypothetically contains relic axions 
born in the Early Universe, and in the cosmological 
context the pseudoscalar (axion) field cj) can be consid¬ 
ered as a function of cosmological time only, The 
time variable t corresponds to the following choice of 
the background space-time metric: 


ds^ = dt^ — a^{t) 


dx^^ + dx"^^ + dx^^ 


(35) 


In this model the gradient four-vector is of the 

form Vi(j) = Uif), where Ui = <5° is the global velocity 
four-vector. In the context of cosmological application 
we consider the axion field potential V{4i^) to be of the 
form where TO(a) is the axion mass. 

Then the equation of the axion field evolution (1551) is 


^ + 3-^ + mf )4> = 0 , (36) 

where the dot denotes a derivative with respect to time. 

We consider the space-time background to be fixed 
by the corresponding gravity field equations, and the 
scale factor a(t) is a known function of time [^ . For 
instance, in the de Sitter-type regime of cosmologi¬ 
cal expansion the scale factor is of the form a{t) = 
a(7o)exp [Ho{t—to)] with the constant Hubble function 
Fl{t) = ^ = iJo- When m(a) > |7Jo the solution to 
Eq. (|551) is 


= e {())(7o)cosH(a)(t-7o)+ 




(a) . 




sinH(a)(t-to) U (37) 
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where the effective axionic frequency is introduced as 


n 


a) — 



(38) 


In these terms the equations for the transverse compo¬ 
nents take the form 

s+ = -n{t)S -, S- = n{t)s+, ( 47 ) 


When the electromagnetic field is absent, the 
equations of particle dynamics (l32l) for the metric (|35]) 
are reduced to 


dPj 

dr 


2m 




(39) 


where 


_ ^257*3(^0) Ht ) 

^ ’ a{t)pOit) ■ 

The solutions to (lT7)l are of harmonic type 


The solution is known to be the following: 


5+= 5 cos (t), S'-= 5 sin'k(t), 


(48) 


(49) 


Pl{t) = Pl{to) , P2{t) = P2{to) , P3{t) = Psito) , 


Po{t) 



(40) 


where = Pi(io)+P 2 (^o)+ 7 ’ 3 (^o) is the constant quan¬ 
tity. Since the space-time is spatially isotropic, we as¬ 
sume that the particle had only one non-vanishing com¬ 
ponent at t = to; say, p^it^) ^ 0 , and consider (for the 
illustration) the following initial data: 

Pi(to)=0, p 2 (io) = 0, S^(to)=0. (41) 


Clearly, for this case at an arbitrary time moment 
Pi{t) = P2{i) = 0 and P 3 {t) = P3{to), so that the cosmo¬ 
logical time t and the proper time r along the particle 
world-line, are linked by the relationship 

- - (42) 

Jto \/rn?a^{t) + 

The equation of spin dynamics (I32p is now of the form 


^ = —(43) 
Ut m 

Its right-hand side does not equal to zero only for two 
values of indices: i=l and i=2. This means that us¬ 
ing the property S’^pk=S'^pQ-\-S^p 3 =Q we can write the 
equations for the components and as follows: 


dS^ 

dt 


HoS^ 1 + 


pI 


= 0 , = 


S^P3 

Pq 


(44) 


For the initial value 5^(to) = 0 these equations have 
the only trivial solution S^(t) = 0, S^{t) = 0. In other 
words the transverse spin components only are influ¬ 
enced by the axion environment and the longitudinal 
component is not touched. 

From the normalization condition S^Sk = —S^ 


we obtain a -I- -S'!) = 5^, which is 

introduce two convenient variables 

the hint to 


a[t) a[t) 

(45) 

so that 

+ S'^ = = const. 

(46) 


where the phase of rotation is presented as a formal 
integral 


^'(t) = [ n{t)dt + 'i/{to) . (50) 

Jto 

For the illustration of the obtained exact solution we as¬ 
sume, first, that the particle is ultrarelativistic (q^ >> 
m?a^{t)) and Psito) is positive; second, we choose the 
time moment to so that 4>(to) = 0; third, we put for 
simplicity 4t(to)=W25()'(to)- Then, keeping in mind that 
nowadays m(a) ^ Hq, we find the phase 'l'(t) in the 
explicit form (see (l37l) I 

^'(t) = Cr cos m(a)(t - to), a = UJ23(t>(to) ■ (51) 

The solutions for the spin components can be now pre¬ 
sented as follows 

— S C 0 S{(TC 0 S [m(a)(t — to)]} = 

OO 

= Jo(cr)-72 ^(-1)" J 2 „(cr) cos [2nm(a)(t-to)], (52) 

n —1 

S- = S sin{tTCOsm(a)(t — to)]} = 

OO 

= 2 ^(-1)" J 2 „+i(cr) COS [(2n-7l)m(a)(t-to)], (53) 

n —0 

where are the Bessel functions of the first kind. 

We deal with sophisticated spin precession, for which 
the phase of precession 'l'(t) is the harmonic function 
oscillating with the axionic frequency n(a)=W(a) and 
the amplitude cr = ui 234 >{to)- 

B. Spin precession of relativistic charged particle 
in storage rings 

When a relativistic charged spin particle moves in 
the constant magnetic field B, the motion is known to 
be circular, and the quantity reciprocal to the Larmor 
frequency, predetermines the time scale of 

dynamic processes. In fact, for such motion the cos¬ 
mological phenomena can be considered as extremely 
slow, and one can put H{t) —>• 0. The scale factor a{t) 
can be replaced by constant a{to) and absorbed into the 
redefined coordinates (in fact, one can put a{t) —>■ 1). 











Since the size of the storage ring is much smaller than 
the typical size of the dark matter inhomogeneity, we 
can neglect a spatial dependence and consider the axion 
field as a function of time only. In our case the quantity 
(j) can be expressed in terms of the energy density scalar 
kF(a) and pressure P(a) attributed to the axionic dark 
matter as follows (see, e.g., [^1: 

^ + -P(a)(0 • (54) 

When the axionic dark matter is cold, i.e., P(a)=0, and 
W(a)=/9(a)) where p(a) is the mass density of the dark 
matter, this formula is simplified, respectively, as 

'^ = ■ (55) 

Also, for the sake of simplicity, we can neglect the defor¬ 
mation of the initial magnetic field by the axionic field, 
and put equal to zero all the new coupling constants 
except a; 25 . 

Let the magnetic field be directed along the = z 
axis, i.e., only one component of the Maxwell tensor, 
Fi 2 = const, is non-vanishing. In the cylindrical coor¬ 
dinates {p, if, z} with the metric 

ds^ = dt^ — {dz^ + dp^ + p^dp^) (56) 


the equation of particle motion 


Dpj 


= -FjkP^ , 


Dt m 
can be transformed as 

dt po pq 

and gives the evident solutions 

Pz{t) = PziO) = 0 , pp{t)=0, 

p^ = -epFp^ = -ep^Fi2 , 

p{t) = R = const, p{t) = , 


Po = \lrrF+e^R'^Fl^ > ^ = ''"l/ ^ 




(57) 


(58) 


(59) 

(60) 
(61) 

(62) 


Here R is the radius of the circular orbit, and the quan¬ 
tity ^(B) given by 


^{B) 


eFi2 

-= const, 

Po 


(63) 


is the relativistic angular frequency of rotation (the rel¬ 
ativistic Larmor frequency). 

The equations for the spin evolution 


^ ^ -F\S>^ (64) 

Dt m m 

rewritten as 

P0^ + — = ^^E^^°^SkPv+eF\S^ 

at p ^ ^ p 

(65) 


give 


= 0, 5“^ = 0 , 

(66) 

= a;250H(B)7? RP , 

(67) 

SP = -uj2o^n(ji)R . 

(68) 

Physically motivated solutions to these equations are 

O 

II 

S- 

o 

II 

o 

(69) 

S^{t) = -5cos4'(H)(t), 

(70) 

SP{t) =5sin4'(H)(t), 

(71) 

where the hybrid precession phase 'L (h) 

is given by 


4'(H)(t)=«'(t) L!(b)P = «'(t)- 


eFi2R 


and '!'(<) is the axionic phase (we put here to = 0) 

4'(t) = a;25[(('(t) - <(>(0)] = 


(72) 


= W25 


0(O)(cosII(a)t-l)+ Sinfl(a)t 


(a) 


~ W25 <^(0) t . 


(73) 


Clearly, the spin four-vector is orthogonal to the parti¬ 
cle momentum four-vector, i.e., pkS'^ = 0. This solution 
can be illustrated as follows. If ^ 25 '/' = 0, the particle 
has the spin three-vector directed along the magnetic 
held, and this direction is conserved during the particle 
circular motion. If uj254’ 7^ 0 fh® spins start to process 
in the plane pOz according the law described by formu¬ 
las dMD-dZIl); the frequency of the precession depends 
on the particle energy, or, equivalently, on the orbit ra¬ 
dius R. From the experimental point of view, if the 
polarized beam of electrons is formed in a storage ring, 
and all the spins are initially directed perpendicularly 
to the ring plane, one can expect, that the axionically 
induced spin rotation will start, and the distribution 
of the angles between the spins and the ring plane will 
be a predicted function of time and the particle energy. 
In the ultrarelativistic regime, when m —>■ 0 effectively, 
one obtains from (17^ that 4'(H)(t) —>■ 'l'(t) and the de¬ 
pendence on the parameter R disappears (some specihc 
details of ultrarelativistic spin particle motion can be 
found also in [^1. 


C. Spin dynamics in the field of an axion star 

In this application we consider a spherically sym¬ 
metric static axionically active object, which is charac¬ 
terized by the metric 

ds^ = B{r)dt^ - A(r)dr^ - {dO^ + sin^ Odp^) . (74) 

The axion field is assumed to depend on the radial coor¬ 
dinate r only. The gradient four-vector S/icj) = Sl(f>'{r) 
(the prime denotes the derivative with respect to r) 
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is now space-like, i.e., Vi(j)V''cl) = < 0. When 

we consider the distribution of the pseudoscalar (ax- 
ion) field in this static model, and the gravity field is 
assumed to be strong, it seems to be reasonable to use 
the extended potential 


The function (j){r) satisfies now the equation 


= A{r)(p 



™(a) + 


(75) 


(76) 


which is (due to the absence of electromagnetic field) 
formally the same as for the gravitating scalar field 
$(r) (see, e.g., [70l - [^ l. There is no need to present 
explicit solutions of (1751) for our purposes. Examples of 
solutions with asymptotically flat space-time and scalar 
fields vanishing at r —^ 00 can be found in [toI - I^ . 

In the static spherically symmetric case the plane 
of the particle motion can be chosen as the equatorial 


one 9 = and the first equation in 
give four integrals of motion 73| : 


is known to 


Po=K=cou.st , pe = 0 , Pip = —J = const, (77) 

1 /IF2 J2 

For the asymptotically flat space-time with A(oo)=l 
and 5(00)=1 the normalization condition pip^=m? 
yields 


K = po{oo) = y'm‘^+Pr‘^{oo ), (79) 

i.e., the constant K usually describes the particle energy 
at infinity. The constant J relates to the conserved 
particle angular momentum. The constant E regulates 
the parameter along the particle world-line, e.g., when 
dr = ds, we see that E = m. 

In order to analyze properly the equations of spin 
evolution (1321) we distinguish two specific types of parti¬ 
cle motion: first, the radial motion; second, the motion 
along a circular orbit. 


5% + S-^Pr = 0, we can write this subset of equations 
as follows: 


dS° 

dr 


+ S' 
dS^ 

dr 


(82) 


2m V A 


B 


As in the first application we assume that 5'’'(0) = 0, 
providing that 5'°(0) = 0 from the orthogonality condi¬ 
tion. Then the solutions to the equations (|8T1) and (|82l) 
are trivial S^{t) = 0 and S'’'(t) = 0. 

The second subset of (l80l) 


dS^ 

dS^ 


dr 


sO ( ^ 1 = - 

mr 


+ — = 


^254>' 
m VABr- 


rPoSip , 


mr 




my/ ABr 


~P0^B ; 


can be transformed into 

using the relation between r and r (see (1781) 1 


T = / dr 


m^J AB{r) 
y/K^ - E^B^ ’ 


and the following definitions: 


5'+ = rS^ , S- = rS'^ , 

- E^B{r) 

Clearly, the solutions to (1551) are 

5*+(r) = 5 cos d/(r), S'-= 5 sin'k(r), 

'l'(r) ='I'(oo) -I- f dr^l{r). 


(83) 

(84) 

(85) 


( 86 ) 

(87) 

( 88 ) 

(89) 

(90) 


Thus, when the particle moves in the radial direction, 
one deals with a spin turn in the plane (0, (j)); the quan¬ 
tity r2(r) plays the role of the rate of turn with respect 
to radial variable r, and 4'(r) describes the cumulative 
angle of turn. 


2. Circular particle motion 


1. Radial particle motion 

For this type of motion pg=0, Pip=0 and the equa¬ 
tions of spin evolution 

^ = ^cj>\r) e*'=’'°SfcPo (80) 

Ut m 

can be split into two independent subsets. The first 
subset contains the components S° and S’’ only, and 
does not include the coupling term proportional to a;25. 
Keeping in mind that the condition S^pk = 0 leads to 


The circular motion is characterized by r = 
R =const, so that p^ — 0 and thus 





(91) 


We consider a stable orbit and thus the radial com¬ 
ponent of the gravitational force should vanish on the 
orbit, providing = 0, or 


B'{R) 2J2 

B^{R) ~ K^R^ ■ 


(92) 



























10 


Also, as previously, we have that pq = K, pg = 0, p^p = 
— J, but now the particle can not reach infinity, and 
we have to re-define the constant K. For instance, we 
obtain from dOTI) and ((Ml) that 




E^R^B'iR) 
2B - RB'{R) ’ 


K = Ex 




2B - RB'{R) 


(93) 


The first equation of (IMl) gives implicitly the radius of 
the circular orbit as a function of orbital moment, R{J). 
The second equation defines the energy of particle on 
the given orbit as a function of obtained radius, K{R). 
We assume that the inequality 2B{R) > RB'{R) is sat¬ 
isfied (for instance, for the Schwarzschild metric B{r) = 
1 — this inequality means that R> = 3GM). 

The orthogonality condition S^pk = 0 and the in¬ 
tegrals of motion dzzl) provide that the components 
and S'^ are proportional one to another, = S'^ {^)- 
For the circular orbit one can replace the differentia¬ 
tion with respect to proper time r with the azimuthal 
angle ip due to the relationship ^ = (recall, 

that in the case s = r we obtain E = m). With this 
replacement three independent equations for spin evo¬ 
lution take the following form: 


The quantity plays the role of frequency of the 

spin turn. Other components of the spin four-vector 
are, respectively: 

[Cl sin Or - C 2 cos Or] + C 3 , 

1[R) 

u9 ( T?\ 

— Tt [^1 ~ ^2 cos Or] -|- C 4 .(100) 

I{R) 

The constants of integration Ci, C 2 , C 3 , C 4 are con¬ 
nected by the normalization condition 

= -5^ = const, ( 101 ) 


which at r = R yields two relationships 


TD Z7’2 p2 

52 ^ (Ci"+C|) +R^cl 




1 -b- 




,PB 


C 3 = -C 4 


iM2^R^(t>'{R)K 

J^AB{R) 


( 102 ) 


On the other hand, three independent constants of in¬ 
tegration Cl, C 2 , C 4 , are connected with initial values 
S'‘^( 0 ), 5®(0), S'’'( 0 ) as follows: 


dS^ 

dp 


fjcS 

S^-H^R), =S^-H\R), 

dip 

jqcp or 

^ + ^ = -S^ -H^iR). (94) 

dp R 


Here, for short, we introduced three auxiliary functions 
of the radius R 


H^{R) 


RE^B{R) 

K‘^A{R) 


H'^{R)=H\R) 


H\R) = 


UJ25R^E^x/Wm^'iR) 

KJx/A(R) 


E‘^B{R) ’ 
(95) 


The evident differential consequence of (IM1) - (IM1) is the 
following equation of the second order for S^'. 

rft QV 

-^+I^{R) = 0, (96) 

d(p^ 

where 


P{R) = 

R 


E^B 

K^A^ J 2 A 


(97) 


For the positive metric functions B{R) > 0 and A(i?) > 
0 the solution to this equation is harmonic function of 
the azimuthal angle 


5'‘^((p) = Cl cos/:/? -I- C 2 sin/tp , (98) 

providing the following solutions in terms of r: 

(t) = Cl cos Hr -I- C 2 sin Hr , 

5'°(r) = [Cl cos Hr -|- C 2 sin Hr] . (99) 

K 


c, = s-(0). C, =-Aa - s«(0). 


RI{R) 

C^ = -S^{0)^+S\t)) 


E‘^IB{R) 


rt 


(103) 


The problem is solved completely. The formulas (IMll - 
(IIOOI) with constants given by (11021) . (11031) and auxil¬ 
iary quantities (Hil), dlTD describe the axionically in¬ 
duced turn of the spin four-vector of the particle mov¬ 
ing along the circular orbit around the static spher¬ 
ically symmetric gravitating object. When a;25:/)'=0, 

we see that I =^\/and the corresponding frequency 


ij 

mR‘^ 


takes the form H = H(geodesic) = 


KR^ 


— de- 
A’ 


H= 

scribing the geodesic precession. This fact allows us to 
indicate the term H(axion) = axionic fre¬ 

quency. With this terminology we can say, that the 
total frequency H satisfies the equality 


H — yf^^geodesic) + ^(axion) ’ (1^^) 

and can be called as a hybrid frequency of the geodesic- 
axionic precession. Finally, it should be mentioned that 
one can introduce local frequency ui instead of H, using 
the equality uidt = fldr. Clearly, we obtain that oj = 

nB{R)f. 


D. Spin precession induced by plane-symmetric 
axion-gravitational waves 

The third application of the model relates to the 
case, when the gradient four-vector is the null one. 
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i.e., \7i(j)V^(j) = 0. It can be realized, e.g. , in the model 
with plane-wave symmetry (^ . [tJ . [^ . The corre¬ 
sponding space-time metric is of the form: 

= 2dudv—L^ [e^^(da:^)^-|-e“^^(da:^)^] , (105) 

where u = and v = are the retarded and 

advanced times, respectively, and two metric functions 
L{u) and /3(u) depend on the retarded time u only. On 
the plane-wave front u = 0 the initial data are fixed in 
the form 


L(0) = 1, L'(0) = 0, /3(0)=0. (106) 


In order to solve the key equations (11091) and (IllOl) . 
it is convenient to use the following auxiliary functions: 


S-{u) = Le^S"^ , S+{u) = Le-^, 

Vl{u) = UJ 2 ^ 4 >'{u) , (113) 


C, 




a 


Ca (e 


a 


In these terms the equations (11091) . (111011 take the form 
^S- = VL{u)S+ -f flg(u) - f'{u ), 

-^S+ = -VL{u)S_ - 0/(u) - g'{u ), (114) 


We assume that the background pseudoscalar (axion) 
field also depends on retarded time only, (j) = (^(u), pro¬ 
viding the condition Vi(j)V^(l) = 0 automatically. Exact 
solutions of this type (in particular, the solution linear 
in the retarded time) can be found in Q. 

The equations of particle dynamics in the metric 
(11051) are known to yield (see, e.g., [^ 1 

= Pv = Cy , P2 =C 2 , P3 =C3, 

m2 + L-2(e-2/5c| + e2/5c|) 

P =Pu = -- ■ (107) 

Here Cy, C 2 , C 3 are constants of integration. 

Our aim is to solve the equations of the spin evo¬ 
lution, which can be reduced now to the following three 
independent equations: 


and their solutions happen to be very simple: 

= A cos 'k (u) — g(u), 

S-= A sin'l'(u) — f(u), (115) 

4'(u) = ^f(O) -h LV 25 WU) - '?^’(0)] • 

Here A is an integration constant. When the integral 
Sy is non-vanishing, the last unknown function Su{u) 
reads 

Su=T^ [f+g‘^- 2 A{gcos'i!+fsin'^)+A‘^-S^] , 

ZDv 

(116) 

and the constant A can be found from the orthogonality 
condition (I112p yielding 


-^5'“ = 0 
du 

e~^ d 


L du 


[Le^S^] 


Ul23(t>' ( 

-iy\ 


d 

L du 


(Le-^S^) + 


W25</>^ 

L2 


= Sy = const 
§-S, = 

(L-Vf)' = 


(108) 


(109) 


( 110 ) 


Here we used the relationship u = t— between the 
parameter r and the retarded time u, which is the con¬ 
sequence of the equation = p^ = Cy (we restrict 
ourselves by the case, when u=0 corresponds to r=0). 
The prime denotes here the derivative with respect to 
retarded time. We present only three equations from 
four, since the component S'"=Su can be found from 
one of the two integrals 


2SuSy = {Le^S^Y + {Le-^S^f - , ( 111 ) 

SuCy + SyPu + S'^C2 + S^C3 = Q. ( 112 ) 


Below we illustrate the obtained exact solutions by the 
examples of longitudinal and transversal particle mo¬ 
tion with respect to the plane front of the gravitational 
wave. 


1. Longitudinal motion 

Let the spinning particle start to move along the 
- axis, i.e., _P2(0) = C 2 = Q, P3(0) = €3 = 0 . Ac¬ 
cording to (11071) at u > 0 the particle keeps the di¬ 
rection of motion, P 2 (u) = P 3 {u) = 0. From (11081) 
and (llll|) it follows that the components Sy and Sy 
do not feel the influence of axions, and we assume that 
^^(O) = S'„(0) = 0, i.e., at u = 0 there were only two 
non-vanishing spin four-vector components, 5^(0) ^ 0 
and 5'^(0) 0. For such initial data we obtain imme¬ 

diately that /(m)= 0, g{u)=0, and thus 

Su{u) = 0 , Sy{u) = 0 , 

-d\ 

— j sin4'(u), (118) 

ad \ 

^ j COS 'k (u). 


S\u) = S\^- 

S^{u) =s(^ 


Clearly, when Si,=0, one can extract S'u from the second Again we deal with axionically induced spin rotation 
relationship only (Cy ^ 0 for massive particles). with the frequency fl{u) = uj 25 (f>'iu)- 
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2. Example of a transversal motion 

Let the particle start to move at it=0 in the direc¬ 
tion (in the front plane of the gravitational wave), 
and have initially only one non-vanishing component of 
the spin four-vector S'^(O) ^ 0. Mathematically this is 
possible if 

S'„=0, «'(0)=0 ^ f(u)=g{u)=Q, A = S. (119) 

Then the exact solution obtained above gives: 

5'„(u)=0, Su{u) =-^ S sin'll(u), 

(u) = S sin'I'(u), (120) 

S'^(u) = 5 cos'L(u). 

The first and second formulas in (11201) yield 

*'("> = (^) ”"*<“> ■ 

( 121 ) 

One can see, that due to the spin-axion coupling the 
longitudinal and the second transversal components of 
the spin appear S^{u > 0) ^ 0, S‘^{u > 0) ^ 0. The 
axionically induced spin rotation is characterized by the 
frequency Q.(u) = a;25(()'(it). 

V. DISCUSSION 

We have established the model of the pseudoscalar 
(axion) field action on the spinning particle. In its 
minimal (curvature independent) version this model in¬ 
cludes the dynamic equation (l30l) and the equation of 
the spin evolution (IHTl) . The terms, which include the 
Maxwell tensor Fmn and its dual are constructed 
phenomenologically by analogy with (and as general¬ 
ization of) the well-known Bargmann-Michel-Telegdi 
(BMT) model. As in the BMT model, the dynamic 
equation (1501) does not contain the spin four-vector, and 
the equation (lOTl) is linear in . Is it possible to ex¬ 
tend m by introducing the spin four-vector quadrat- 
ically, e.g., as it was made by Bander and Yee in 
For sure, the scheme, based on the representations (ED 
of the master equations and on the decomposition of 
the basic tensor introduced in m, gives us such 
tool. For instance, when w** is linear in the spin, the 
corresponding dynamic equation is also linear, and the 
equation of the spin evolution is quadratic in the spin 
four-vector. 

However, even in the simplest BMT-like form (15(11) 
the dynamic equation includes two new force terms, 
which points to the axionic extension of the theory. The 
first novelty is the term in which the tensor eF^n is re¬ 
placed with the tensor eejjFf^^, where the pseudoscalar 


multiplier cj) compensates the pseudo-tensorial nature of 
the dual Maxwell tensor The second novelty is the 
term with the gradient-type multiplier m front of 

. If the axionic environment indeed produces forces 
of these kinds, they could be tested in high-energy ex¬ 
periments with polarized beams, e.g. in the LHC. We 
hope to present the corresponding work and discuss ex¬ 
act solutions to the whole system of equations of axion 
electrodynamics, particle dynamics and spin evolution 
in the next paper. 

As for this paper, we consider the contribution 
from the only new term, which is free of the Maxwell 
tensor and linear in the gradient four-vector This 
term describes the direct action of the axion field on the 
particle spin. This term, SkPm was intro¬ 

duced phenomenologically, and the dimensionless cou¬ 
pling constant 0125 should be recognized. One of the 
ways to find 0125 is to make the reduction from the ax¬ 
ionically extended Dirac theory; we will return to this 
problem in the future. The second way is to use the 
analogy with the axion-photon coupling, which has been 
considered in Section m If to follow the hypothesis of 
universality and to compare the evolutionary equations 
for the photon polarization ED and for the spin rota¬ 
tion ED, one can assume that W25 = a = 1. Neverthe¬ 
less, one should repeat that this coupling constant has 
to be found experimentally. 

One can mention that the term 
describing the direct spin-axion coupling, can be repre¬ 
sented in the BMT-like form. Indeed, let us take the 
main term ■^F^'^Sk appeared in the BMT model for 
the case g = 2, and consider the decomposition of the 
Maxwell tensor in the reference frame associated with 
the particle moving with the velocity In addition, 
let us assume that in this frame the electric field is 
absent, then we obtain that F** = where 

is the four-vector of the corresponding magnetic ex¬ 
citation. Comparing the terms and 

— :^f^^’‘"'SkBiPm we can see that they formally coin¬ 
cide, when 0025^1(1) = —^Bi. In other words, this anal¬ 
ogy hints that the gradient of the pseudoscalar (axion) 
field can produce the spin rotation similar to the well- 
known effect induced by the magnetic field. 

Of course, the mentioned analogy is incomplete, 
since B^ is the space-like four-(pseudo)vector, while 
the gradient four-(pseudo)vector Vfcc/i can be time¬ 
like > 0), space-like {Wcj) < 0) or null 

{WicjiV^(j) = 0). However, in all three cases, as it was 
demonstrated using the obtained exact solutions to the 
master equations, we deal with the same phenomenon, 
the axionically induced spin precession. 

The typical example for the time-like gradient 
four-vector is given by the spatially homogeneous 

cosmological model according to which the pseudoscalar 
field corresponds to the relic dark matter axions. In this 
case the gradient four-vector reduces to S°ip and the 
spin rotates in the plane orthogonal to the direction of 
the particle motion. The corresponding time-dependent 
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axion field (0 ^ 0, 0'(r) ^ 0, 0'(u) ^ 0), activates 
the spin-axion coupling. In this sense, when the grav¬ 
ity field is strong, it displays the phenomenon of spin 


frequency fl(t)=a;250(t)^ (V is the modulus of the ve¬ 
locity three-vector) is a direct analog of the Larmor fre¬ 
quency. 

The static spherically symmetric model of gravi¬ 
tational and axion fields gives the typical example for 
the space-like gradient, Vi0V*0 < 0. When the particle 
moves in radial direction, we deal again with the spin 
rotation in the transverse plane. Now it is more reason¬ 
able to speak about the spin turn with respect to radial 
variable r rather than with respect to time (see (I5S1) - 
(IMll i. The rate of spin turn depends on the distance 
to the center; in particular, when pr{oo) = 0 and thus 


K=m=E, we obtain r2(r)=a;25- 


4''{r 


I , Far from the 

center the metric function B(r) has the standard be¬ 
havior, B{r) —(l— , so the asymptotic behavior 

of the frequency V,{r —>• oo) — {r) is 

predetermined by the function One can expect 

that the zone of strong gravitation gives the maximal 
contribution into the total turn of the spin, however, 
this question should be analyzed in its own right. 

When the particle moves along circular orbit 
around the axionically active object, one can split the 
total effect in the spin rotation into geodesic preces¬ 
sion and axionic precession. The frequency of rotation 
is given by hybrid formula (110411 . and the axionic fre¬ 
quency (axion) = Constant on the orbit, but 

depends on the radius of the orbit R. The behavior of 
the function fI(axion)(-R) can be studied only when the 
distribution of axion field 0(r) is found. We hope to 
return to this question in the next paper in the context 
of discussion of qualitative and numerical study of the 
total system of master equations. 

The case Vi0V*0=O is typical for the model with a 
plane-wave symmetry, for which the axion field depends 
on the retarded time only, 0(u). Again the exact solu¬ 
tions to the master equations demonstrate the spin ro¬ 
tation with the frequency n{u) = W 250 '(m). Two cases 
have to be distinguished in this model. First, when the 
particle moves orthogonally to the front of the gravita¬ 
tional wave (the so-called longitudinal motion), we deal 
with simple spin rotation in the front plane ^ 0 
and 5'^ ^ 0). When the projection of the particle mo¬ 
mentum on the front plane is non-vanishing {p 2 ^ 0, 
transversal motion), the spin rotation becomes more 
sophisticated (the additional, third component of the 
spin four-vector appears). 

The first obvious conclusion for all three applica¬ 
tions, is that the pseudoscalar (axion) field makes the 
space-time chiral, so that the left-hand and right-hand 
rotations of the particle spin four-vector become non¬ 
equivalent. The spin precession can be indicated as the 
first general property of the model. 

The second general property is that the gravita¬ 
tional field, providing the non-vanishing gradient of the 


rotation more effectively. 

Since the phase of the spin turn is described by the 
integral formulas of the type 4' = / {^=t, ^=r 

or ^=u), the effect of spin rotation is cumulative in the 
space-time domains, in which the quantities 0(t), 0^(r) 
or 4>'{u), respectively, hold the sign. In this sense the 
phase accumulation can be treated as the third general 
property of the model. 

Finally, we would like to say a couple of words 
about estimation of the described effect. We prefer 
to do it on the example of relic dark matter axions, 
which seem to be distributed everywhere, using the 
model of relativistic charged spin particle motion in a 
storage ring with magnetic field (see Section IIVB|) . In 
the cosmological context the quantity 0(to) can be es¬ 
timated as 4>{to)=-^\/P(a)ito) using the mass density 
of dark matter axions P(a)(to)- Thus, for the relic cold 
dark matter axions with the mass density of the or¬ 
der P(dm) — 0.033 M(gun)pc“^, for the ultrarelativis- 
tic particle with F —>■ c, for the coupling constant 
■^=PA77 — 10“®GeV“^, we obtain that an optimistic 
estimation for the spin rotation frequency (in Hz) is 

f^(H) = ^(H) ^ ^ 


~ 10-6 ( PAj'r \ f Vf^(DM) 

V 1 V10-6GeVV W1.25GeV- cm-3 

In order to estimate the possible total axionically in¬ 
duced phase variation, Adt, for other examples, we have 
to know the time period, during which the sign of the 
frequency is non-changed and is, say, positive. When 
we deal with homogeneous cosmological model, accord¬ 
ing to (15^ . (1551) this time period is about of T = 
and thus is negligibly small. The application to the 
static spherically symmetric gravitation field is much 
more promising, since now the derivative (j)' (r) is mono¬ 
tonic function, and the accumulation of the phase of the 
spin turn can continue during long time for both radial 
and circular motion of the test particle. 
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